Continuous variable noise-free states in correlated quantum noisy channels 
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We explicitly compute the evolution of the density operator of a two-mode electromagnetic field 
when the influence of the thermal fluctuation of the vacuum is common for both modes. From this 
result, we give an example in which the bundle of quantum noisy channels turns out to be noiseless 
for the special type of signal states due to the correlation. 



Computing the classical capacity of quantum channels is one of the most important issues in quantum information 
theory 1]. The quantum channels are represented by a stochastic or completely positive map that maps the density 
operator tt to the density operator cr; tt i— > cr = $(7r) ^. If the channels $i and <i>2 have no correlation, the bundle 
of them is simply given by $1 (8" *&2- Even in this simple case, whether the classical capacity of quantum channels is 
additive or not is still a long-standing open question. On the other hand, the correlation or memory effect between 
channels leads to interesting phenomena. One such example was given in 3], where the entangled signals enhance the 
mutual information of the pair of depolarizing channels with memory. Without correlation, the depolarizing channels 
are proven to be additive [j| and explicit calculations have shown that the input of entangled qubits cannot increase 
the classical capacity of depolarizing channels 's). Similar results for channels with memory have been obtained in 

Qi 1^ ' Asymptotic classical capacity of quantum channels with memory has been also analyzed Q . 

In this paper, we give an example in which the bundle of quantum noisy channels turns out to be noiseless for the 
special type of signal states due to the correlation. This constitutes one example of the noiseless codes ^3 ^| in a 
continuous variable system. 

Consider a two-mode electromagnetic field with the same frequency luq interacting with the reservoir: the thermal 
fluctuation of the vacuum. If these two modes are spatially well separated, the reservoir field acts independently on 
each mode. The interaction Hamiltonian takes the following form. 

i=l j 

Here, denotes the creation operator of the electromagnetic field of mode i and fej*-* denotes the annihilation operator 
of the electromagnetic field of reservoir i acting on mode i. They obey the commutation relation 

^fc^'^z"'''^ ~ ^ij^ki- Other commutators are zero, gj*'' is the coupling constant. This independent reservoir model 

has been well investigated [l^ 113. Iisj . However, if the two modes are close enough spatially, the action of noise on 
each mode must be more or less correlated. More precisely, the ith reservoir field 6*^*-' acts partially on mode i). 



= Sij and 



— 5ij. Here, bj is the annihilation 



Therefore, in Eq. should be replaced by y/T^^bj -\- y/JIbj with 

operator of the common reservoir that is independent of reservoirs 1 and 2 and the parameter /i (0 < A* < 1) represents 
the degree of correlation. Under the extreme condition of /i — > 1, the interaction Hamiltonian takes the form 

^ = ^(E«I)E-9A + ^-C-, (2) 

where gj = g^^"^ — gf^- Namely, the reservoir acts on both modes simultaneously. The same model has been applied 
to the two-mode squeezed vacuum states . This ideal model may be hard to achieve experimentally, but it provides 
the essential feature of these correlated channels of intermediate value of /i (0 < /i < 1). Hereafter, we confine ourselves 
to this extreme case. 

Following the standard procedure for deriving the master equation for the density operator under the Markov 
approximation in the weak coupling limit |l7|. we get the following equation of motion for the two-mode density 



2 



operator in the interaction picture. 



pit) = 77(^0 + 1) 2{ai+a2)pit){a\+al)-{a\+al){ai+a2)pit)~ p{t){a\+al){ai+a2) 



dt 



r 



2{a\ + a3,)p(t)(ai + a2) — (ai + a2){a\ + a|)p(t) — p{t)(ai + a2)(a| + aj) 



(3) 



where iVo = 1/ {e^"^"/^^ - l) is the mean number of reservoir quanta at frequency ojq. The initial [t = 0) density 
operator can be expanded in terms of coherent states as follows. 



P(0) 



■ j W d^ad^l3,P (ai, as; PI, (3*2) A ("i: "2; PI, 0*2) , 



(4) 



where 



2 

y"--- J f[d^a^d^(3,P {ai,a2;f3l,f};) ^ I, 



(5) 



and 



A(ai,a2;/3r,/32*) 



I«l,a2) (/3l,/32| 

(/3i|ai)(/32|a2) 



J|exp (-tti/?*) exp(Q;ja|) |0) (0| exp(/3*ai) 



(6) 



is an analytic function of ai, a2, Pi, and P2- p(t) is the output state for the input state p(0) in the channel with 
correlation considered here, $c; p{t) = <I>c [p(0)]- 

In order to solve Eq. (|3J) under the initial condition Q), we first calculate the two- mode characteristic function 
defined as 



JJ^exp (Xia] - A+flij p{t) 



x(Ai, A2;i) = Tr 

Following the standard procedure 0|, we get the equation of motion for x(Ai, A2; from Eqs. (O and (0) as 



df 



x(Ai,A2;t) = -rnVo + - |A+rx(Ai,A2;0-r a;— + A+— x(Ai,A2;t), 



d 



d 



dX^ 



where A± = Ai ± A2. From Eqs. Q and (jT)), the initial value of x(Ai, A2;i) is calculated as 



x(Ai,A2;0) = 



Y[ d^ad^P.P {ai,a2;Pl,Pl)W_ exp 



X*ai + XiP* 



(7) 



(8) 



(9) 



Following the operator method described in [l^, we can solve Eq. © under the initial condition ©. The straight- 
forward calculations yield 



x(Ai, A2;t) = 



\{d^ad''P^P{al,a2■,Pl,P*2) 



X exp 



-\N{t) |A+|^ - i |A+|^ - i |A_|^ - ie-r* Ka; - /3;A+) - \ (a_A*_ - P*_X^ 



(10) 



with a± — ai ± a2, /?± = /3i ± /32, and N{t) — Nq{\ — e^^'"*). It is an easy task to check that x(Ai, A2;i) obtained 
above satisfies Eqs. © and ©. 

In the following, we calculate the two-mode Q-function from x(Ai,A2;t) and then calculate the density operator 
from the two-mode Q-function, which is defined as the diagonal element of the density operator in a coherent state; 
Q((5i,(52;i) = {5i,52\ p{t) \5i,52) /tt^. This is also calculated from the two-mode anti-normally ordered characteristic 
function 



X^(Ai,A2;t)-Tr 



exp(-A*ai) exp(Aia]')/9(t) 



(11) 
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via the following integral transformation 

2 

^ •' i=i 

Since x^{^i^^2\t) = e"l'^il''/^"l^2l^/^x(Ai, A2; i), we can calculate the Q-function by using the identity [l7 l ll9 t 



— j d^T^exp A |?7|^ + /ij] + 



V7] J = — CXp 



T 



which holds for Rc(A) > and arbitrary /i and v. The result is 

2 



1 



1 



7r2 1 + N{t) 



(12) 



(13) 



(14) 



where 



f {ai,a2; PI, f3l] 51,82) 
'^ + N{t) f,^ ,2 , ,2 



1 + N{t) 



+ a- + 



Njt) 
1 + N{t) 



{6,5*2 + StS2) - 



1 



1 + N{t) 



Pie 



1 + N{t) 



+ P-]5i + 



1 + N{t) 

Now, we get the density operator from the following formula. 

2 n 2 



1 + Nit) 



- a- U; + 



pie- 



1 + Nit) 



Pit)^[^) /•••/n^^M%exp(-i 



-- |<5, - 7,r ) A (5i, (52; 71*, 72*) Q { ^{Si + li) , ^{62 + I2) it 



This is an obvious extension of the relation 12 C 



1 

P= - 



d^-fd'^SAi-f + S,-f-6)e~^^\ g(7) 



(15) 



(16) 



(17) 



to the two-mode case. In Eq. ((TTIl . A(q!, P) = |a) (/3*| / (/3*| a) and (5(7) — (7I p I7) /tt is the single-mode Q-function. 
At zero temperature, i.e., iVo = 0, the integral of Eq. H16|l . can be explicitly performed by using the identity 



i J d^aexp (- \a\^ + Pa*) /(a) = fiP) 



for an arbitrary analytic function /. The result is 
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pit)^ J ■■■ Jf[d^a,d^P,Piai,a2;Pl,P*2)Aia+it),a.it);P*+it),P*_it)), 



(18) 



(19) 



where a±(i) = (a+e^'"* ± a-)/2 and /3±(i) = (/S+e^'"* ± /3_)/2. It is quite easy to see that p(t) thus obtained 
satisfies the master equation (jJl with iVo = 0. It should be noted that the initial pure state p(0) = |q;i,q;2) (ai,Q;2| 
remains pure; it becomes pit) = <I>c [p(0)] = |a+(t), (i)) (Q;+(t), Q;_(i)| as the output of the channel $c. More 
interestingly, the state |a, —a) (a, —a\ is not affected by the noise of the channel $c. This is also true for the entangled 
coherent state [^ll^ of the form —a) ± e"^ |— a, a)) /\/2. These perfectly decoherence-free states arise from the 
quantum interference coming from the indistinguishability of modes 1 and 2 for the common reservoir. For the input 
state piO) — \a,—a) {a,—a\, the quantum interference is completely constructive while for p(0) = \o:,a) (a, a|, it is 
destructive. 

At finite temperatures T > 0, the input pure state does not keep its purity during the time evolution: the output 
state is mixed even if the input state is pure. The output purity Trp^(i) is calculated from Eq. (|16|l with Eqs. H14() 
and H15|l by using identity (|13|l . The explicit calculations yield 



Trp^t) 



2(1 + A^(0) 
2 + 6Nit) + 57V(t)2 



(20) 
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This is a monotonically decreasing function of t. Another measure for characterizing the channel is the fideUty that 
quantifies how close the output state p{t) and the input state p(0) = |q;i,q;2) {(Xi,ct2\ are: 

F(ai,a2) =Tr[p(OV(0]. (21) 

This is also calculated from Eq. llt)|) with Eqs. H14|) and (|15|l by using identity H13|l . The fidelity for the initial pure 
state p(0) = |a, —a) (a, —a\ is 

F(a,-a)^ 5-—, (22) 



that is independent of a and tends to 1 for A^o ~^ 0. This clearly contrasts with the fidelity 

(23) 



1 + N{t) 



calculated for the initial state p(0) — \a,a) {a,a\. The state |q!, — a) (a, — aj is almost noise-free even for finite 
temperatures while the state |a, a) {a, a\ is much more influenced by the noise. 

In summary, we gave an example of the noise-free or decoherence-free states in quantum noisy channels with 
correlation by analyzing two electromagnetic fields coupled to a common thermal reservoir. The results suggest that 
we can construct a partially noise-free signal state in the bundle of electromagnetic field channels if the modes of 
signal field are spatially close enough to share partially the common thermal reservoir. 
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